On the nonequilibrium relation between potential and stationary distribution 
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We investigate the relation between an applied potential and the corresponding stationary state 
occupation for nonequilibrium and overdamped diffusion processes. This relation typically becomes 
long ranged resulting in global changes for the relative density when the potential is locally per- 
turbed, and inversely, we find that the potential needs to be wholly rearranged for the purpose of 
creating a locally changed density. The direct question, determining the density as a function of the 
potential, comes under the response theory out of equilibrium. The inverse problem of determining 
the potential that produces a given stationary distribution naturally arises in the study of dynam- 
ical fluctuations. This link to the fluctuation theory results in a variational characterization of the 
stationary density upon a given potential and vice versa. 
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I. INTRODUCTION 



Imagine independent colloidal particles in a potential 
field and subject to friction and noise as imposed by a 
thermal reservoir or background fluid. In thermal equilib- 
rium at inverse temperature (3, Prob[a;] oc exp{— (3V (x)) 
where V is the potential on the states x. Typical ex- 
amples include the Laplace barometric formula but also 
the distribution of particles in a fluid undergoing rigid 
rotation. We now add an external forcing and we wait 
till a steady regime gets installed. The stationary statis- 
tics depends on the potential, but most certainly and 
because of the forcing the resulting time-invariant dis- 
tribution of velocities and positions of the particles gets 
modified with respect to the Maxwell-Boltzmann statis- 
tics. The relation between potential and stationary dis- 
tribution is far from understood for generic nonequilib- 
rium systems, beyond its general specification as being 
for example a solution to the time-independent Fokker- 
Planck-Smoluchowski equation. 

We show here that small local variations of the po- 
tential can globally affect the relative density, provided 
a nongradient driving is present. That effect already 
arises in linear order around equilibrium. For the inverse 
relation, we are asking to reconstruct the potential 
which, under a known driving, realizes a stationary 
distribution. In equilibrium, the change in the potential 
needed to change the density p — > p exp(— ^) is simply 
equal to that A. Not so in nonequilibrium! It is then 
a genuine problem what potential field can produce a 
given particle density at a fixed driving. Beyond obvious 
applications in interpreting observational data, this 
question also naturally emerges in dynamical fluctuation 
theory, cf. [2lj. In particular, the large fluctuations of 
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certain time averages around their stationary values are 
governed by a functional (sometimes called an effective 
potential) that is given in terms of the potential of the 
inverse problem. We will give two possible approaches 
to finding the potential, one of which is analytical and 
the other one is based on a variational formula. 



The specific examples we treat in this paper and for 
which the analysis can be applied are those of diffusing 
particles in a background rotational velocity field. Dif- 
fusion in rotating media is one of the central objects in 
geophysical and astrophysical applications. The ques- 
tion of nonlocal and irreversible effects is of particular 
interest for galactic dynamics, where according to Chan- 
drasekhar's theory the huge number of relatively small- 
size gravitational encounters gives rise to an effective 
Brownian motion for test stars, cf. see also [1] for an 
account in the context of galaxy formation. Arguably, 
also experimentally a most realistic scenario for main- 
taining a constant nonconservative force is via differential 
rotation. One can think of concentric cylinders rotating 
at different angular frequencies which are imposed on the 
fluid by stick boundary conditions. Far from equilibrium 
and under a general angular (possibly angle-dependent) 
driving one can expect not only that currents are being 
maintained but also that the time-symmetric aspects of 
diffusion can be essentially changed. It is known that the 
diffusion phenomenon itself may be influenced by rota- 
tion, even by rigid rotation [30| . Related considerations 
on nonequilibrium diffusion also apply to other scenar- 
ios, including e.g. shear flow Ull, M, [m and oscilla- 
tory flows [16], and stochastic models of particle trans- 
port in turbulent media necessarily include discussions of 
driven diffusion phenomena, [l^. Colloidal particles in 
harmonic wells and driven by shear flow have been explic- 
itly treated for the violation of the fluctuation-dissipation 
relation in psj . 

The relation between stationary density and poten- 
tial is generically nonlocal or long range but as our 
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particles are mutually independent, that obviously says 
nothing about the generic long range correlations in 
nonequilibrium systems. As we will see however, the 
mathematical mechanisms are not unrelated. Various of 
these nonlocal effects and generic long range correlations 
under nonequilibrium conditions have been discussed in 
similar contexts of interacting particle systems, most 
recently from the point of fluctuation theory [U, 0, Q , in 
perturbation theory for the stationary density [l^, and 
as a result of the breaking of the fluctuation-dissipation 
relation [l^, [H, [l^l . It was originally the mode-coupling 
theory in hydrodynamic studies for a fluid not in thermal 
equilibrium that revealed the (macroscopic) long-range 
correlations between fluctuations. Light scattering 
experiments can reveal these correlations, [sl ITtI. [26l. [27l| . 
In contrast, the present work analyzes the nonlocality on 
a mesoscopic length scale and that does not result from 
particle interactions but as is already present in a single 
particle distribution (or, equivalently, in a density or an 
ensemble of independent particles), due to the imposed 
nonequilibrium driving. 

In the next Section [ll] we specify our dynamical model: 
independent particles undergoing an overdamped diffu- 
sive motion in a conflning potential and under driving. 
The nonlocality in the relation between potential and 
distribution is discussed throughout Section UIIl The lat- 
ter already speaks about the inverse problem of deter- 
mining the potential for a given density, which is then 
elaborated on in Section [TVl in the context of variational 
principles and for the purpose of dynamical fluctuation 
theory. In the Appendix we add further details about 
the McLennan's interpretation of nonequilibrium distri- 
butions, about the Green's functions encountered in the 
response problem, and about the dynamical fluctuation 
origin of the variational principles under consideration. 

II. DIFFUSION IN A TWO-DIMENSIONAL 
ROTATIONAL FLUID 

Restricting ourselves for the moment to the two- 
dimensional plane with points labeled by the polar co- 
ordinates X — (r, 9) , we consider an ensemble of indepen- 
dent test particles subject to a rotation-symmetric force 
with potential U{r), sufflciently confining so that 

r+oo 

Z^2tt e~^rdr<+(X) 
Jo 

The particles are suspended in a nonequilibrium fluid ex- 
erting an additional force that can have some conserva- 
tive component with potential <I>(r, 9) and a nonconserva- 
tive force v — {vr, vg) for which we assume that the radial 
component Vr = vanishes. The angular driving force 
Vg can be associated with the local velocity of the back- 
ground fluid which is maintained in a differential rotation 
state. We assume a thermally homogeneous background 
(setting the temperature to one) , modeled by a Gaussian 



and temporally white noise. Given that the motion is 
nonincrtial, it satisfies the Langevin equation 

dxt = udt- V(C/-H$)di + \/2dBf (II.l) 

with Bt standard two-dimensional Brownian motion. 

An implicit assumption in the above construction is 
the smoothness of the functions U, <&, and vg on the 
domain taken here as the entire two-dimensional plane. 
Yet, interesting modifications arise when the origin r = 
is not accessible and the particles can only move in 
the non-simply connected domain obtained by removing 
the latter. This allows for possible singularities when 
approaching the origin and the existence of a potential for 
the driving force v is no longer equivalent to the condition 
V X t) = 0; the rotational field of the form v = {0,vg), 
Vg (X l/r serving as example. The exclusion of the origin 
can be ensured e.g. by an infinitely repellent potential 
therein or via suitable boundary conditions at the origin. 

Note that in the absence of a potential, U — ^ — 0, the 
dynamics pi.l|) does continue to make sense, although 
a normalizable stationary distribution no longer exists. 
The transient regime is still relevant and has been studied 
in detail [ssj . 

The stationary distribution for the dynamics (jll.l[) 
has density p verifying the Fokker-Planck-Smoluchowski 
equation 

V-J = 0, J = pv - p\7{U + ^) -Vp {11.2) 

We refer also to [l^, |3l| for more thermodynamic and ki- 
netic gas considerations in the derivation of that nonequi- 
librium dynamical equation. In polar coordinates the 
probability current J = (Jr, Jg) takes the form 

a(C/-f$) dp ^ pd<P Idp 

'Jr = -P ^ -TT, Jg^pvg — — 

or or r o9 r o9 

(n.3) 

By turning on the driving vg, typically not only an an- 
gular current Jg is generated but also a nonzero radial 
component Jr does get maintained. This is a priori not 
in contradiction with the existence of a stationary dis- 
tribution; its normalizability essentially depends on the 
imposed potential [/ + $ and on the boundary conditions. 

Our general aim is to analyze the relation between test 
potentials $ and stationary densities p, under a given 
confining rotation-symmetric potential U and as medi- 
ated by the rotational field vg. First, we examine the 
issue of spatial nonlocality. 

III. LONG-RANGE RESPONSE TO CHANGING 
THE POTENTIAL 

Nonlocal features have been widely discussed in 
the nonequilibrium literature. Mostly however that 
deals with the presence of long range correlations, 
cf. 0, d, 13, m, [131 for time-separated viewpoints, 
or with models of self-organized criticality, cf. pj|. In 
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our case, we have independent particles, hence there are 
no correlations between the particles and correlations 
between spatial points only appear because of fixing the 
number of particles or by fixing the mass. We think of 
the spatial dependence in the density as it is affected 
by local changes in the external potential, and vice versa. 

In the absence of driving, v — 0, the stationarity equa- 
tion pi.l|) has the usual equilibrium solution p cx e^'^^*, 
J = 0, which is manifestly a local functional of the test 
potentials $ in the sense that the response 



log 



p{x) 



~ S{z — y) — S{z — x) 



is insensitive to perturbing $ away from both points x 
and y. 

One can ask to what extent this is an equilibrium prop- 
erty but, in fact, it is easy to devise special nonequilib- 
rium conditions where such a locality still holds. As a 
simple example, take <& = and let the angular driving 
velocity be rotationally symmetric, vg — u{r). Then, the 
stationary density and current become 



1 



J 



(IILl) 



Although the stationary density coincides with the equi- 
librium solution and is perfectly local in the potential [/, 
it now corresponds to a current-carrying steady state. 
Thus one cannot unambiguously decide just from the 
stationary distribution itself whether the system rests in 
equilibrium or whether irreversible flows are present. 

Yet, generic nonequilibrium distributions do get modi- 
fied due to driving and, as a result, they typically pick up 
some nonlocality. Next comes a simple demonstration. 



A. An exactly solvable model 

In the example above the radial currents are absent 
and the steady state is rotation-symmetric. A simple ex- 
actly solvable case where the rotation symmetry is bro- 
ken can be obtained for an angular driving of the form 
ve = f{S)/r and for test potentials that are constant 
along radials, $(r, 0) = ^{6). The origin is excluded 
by the boundary condition Jr{0) — 0. In this case, the 
stationary density is found to be of the form 



p{r,d) ^ p{r)q{9), p(r) 



1 



-Uir) 



which is under the given assumptions the general form for 
a density with everywhere vanishing radial current, Jj. = 
0. The steady state decomposes into separated concen- 
tric motions and the angular distribution q{9) in pil.2p is 
determined from the stationarity condition (jll.2p . which 
reads Jg{r,9) = j{r) for some rotation-symmetric func- 
tion j{r) that can be determined. Explicitly, from pi.3p . 



{f{e)-^'i0)}q{e)-q'i0) 



rjir 



p{r) 



(III.3) 



which decouples the polar coordinates and confirms the 
Ansatz pil.2|) . The angular current j(r) is obtained by 
dividing pll.3|l by q and integrating over the angle vari- 
able, which yields, always for ve{r,9) = f{9)/r, 



P{r) lo f{9)d9 
r f^\-H9)d9 



(III.4) 



As expected, a nonzero steady current is maintained 
whenever the angular driving does not allow for a po- 
tential, i.e., if the work performed over concentric circles 
is nonzero, vu = J^^ f d9 ^ 0. Using the normalization 

condition J^^ q{9)d9 ~ 1, the solution of pil.3p is ob- 
tained in the form 



q{9)^^j\'^^''-'^d9\ 



(III.5) 



n = 



2tt r^TT 



with the work function 



Wie',e) ^^{9')-^{9) 



Here we have used the notation ^g, for the integral per- 
formed along the positively oriented path 0' — > 6* on the 
circle, i.e., it coincides with Jg, for 9' < 9 whereas it 
equals to Jg^ + Jq otherwise. In the sequel we also em- 
ploy the shorthand § for the integral J^^ . 

Equivalently, the same solution to pil.3[) can also be 
obtained in terms of the current j(r)] this leads to the 
next explicit expression for the latter: 



(III.6) 

The equilibrium angular distribution is recovered for 
w = (or j — 0), in which case the work function 
W{9' , 9) derives from a potential and formula piI.Sp boils 
down to the Boltzmann-Gibbs form. 

For w 7^ the character of the stationary density be- 
comes modified, as can be read from the response to 
changes in the test potential $ = '^{9) that we take 
to depend on the angle 9 only. We take the functional 
derivative for changes in the value of ^' at fixed angle rj: 



log 



p{r,9) 



Yir],9')-Y{ri,9) (III.7) 



(III. 2) where 



Y{r^,9)^Sifj-9) 



It is the second term that generates the nonlocality as its 
reciprocal 



^*(r,)-*(',,') + f,;'' {f-w S(i-9)}di 
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still depends on 9. Note also that since the work function 
W{-q,9) is discontinuous at ^ = 77, so is the nonlocal 
contribution in the response (jIII.7p a± -q — and -q — 6' . 

That nonlocal term is manifestly a correction of or- 
der 0{'w), hence, some more explicit information can be 
obtained within the weak driving approximation, consid- 
ering the driving force f{9) (or total work per cycle w) 
small. This yields, up to 0{w): 



= Si-q - 9') - 5{r^ - ^ 




, (7* d^ if ?/ e 
-g* d^ if r/ e 
where (77) stands for the auxiliary density 



and 776 (0 — > 0') indicates that r\ belongs to the posi- 
tively oriented path from 9 to 9' . This explicitly confirms 
the nonlocal and discontinuous structure of the response: 
The relative density p{r,9) / p{r,9') of our weakly driven 
system remains (strongly) sensitive to locally modifying 
the potential, ^{9) ^'(0) + eS(9 — 77), at an arbitrary 
angle rj. The effect disappears at equilibrium, w — 0. 

Note that this special example is essentially one- 
dimensional and we have so far only considered the re- 
sponse to an r— constant change in the potential. The 
response to a strictly local perturbation will be discussed 
in the next section by a more general method. 

The nonlocality of the inverse problem for this specific 
example, potential as function of the density, is continued 
around equations (IIII.13j) - (IIII.15j) . 



B. General linear theory 

Having observed that the nonlocal features of nonequi- 
librium states already emerge in the lowest order of the 
driving strength, we can now follow the linear analysis 
more systematically. The method is by now well known; 
we refer to earlier work of McLennan and others, [24, 36*1. 
The analysis starts from the overdamped form pi.ip , and 
again for simplicity we keep in mind diffusion in the en- 
tire plane under natural boundary conditions (decay at 
infinity) . 

One takes the reference equilibrium distribution (for 
V — 0) as 



Po(x) = ie-^(-) 



(in.9) 



where we combine V = U +^ and Z is the normalization. 
Assume now that the driving velocity field v is uniformly 
small, \v\ = 0{s). 

The solution of the stationarity equation (|IL2p to lin- 
ear order in e reads 

log — = — (V • w - V • VV) (IILIO) 

Po Lv 



where Ly = — Vl^ • V -I- A, which is recognized as the 
generator of a reversible diffusion in potential V\ see Ap- 
pendix |A] for more details and for a physical interpreta- 
tion in terms of the McLennan's theory, [23]. Writing 
equation (jIIL10[) in the form 

Aiy- VF- Vz/ = V-u-w Vy (IILll) 

with V = \og{p/ Po), its variation along V — > V + 5V is 



In terms of the equilibrium generator Ly and the sta- 
tionary current J — p{v — W) — Vp, this reads 



Lv 5v = ---\/5V 



(in. 12) 



always up to terms O(e^); in the same order the density 
p on the right-hand side can be replaced by Po. This 
is a Poisson equation for 5v in which the Laplacian 
(free diffusion) is modified with a drift term in potential 
V . Note that since Lv ^ j-V ■ (poV) and V • J = 0, 
this problem is equivalent to studying the electrostatic 
potential in an inhomogeneous dielectric environment as 
generated by a source with zero total charge. Whereas 
the local component of the response to 5V is already 
hidden in the y— dependence of the equilibrium density 
Po, the nonlocal character of solutions to (jIII.12[) follows 
from general features of elliptic operators. The solution 
of the free Poisson equation (with only the Laplacian 
and natural boundary conditions at infinity) is explicit 
and manifestly long range. The potential V or the 
finiteness of the system introduces an extra confinement 
and the claim needs to be refined. One expects that 
within the confinement region where the density po is 
approximately constant and near its maximal value, the 
response 5v{x) at a; to a local perturbation 6V (both 
localized within that region) can be well estimated by 
replacing Ly with the free Laplacian = A. This 
intuition is indeed correct as we shortly explain in 
Appendix [B] In that region, the response 5v{x) to a local 
perturbation 5V derives from the Green's function for 
d^ dimensional Brownian motion, weighted by the local 
mean velocity J / p. In this sense the nonlocal response is 
intrinsically a nonequilibrium feature. 

Although the above linear analysis can formally be ex- 
tended far from equilibrium and one obtains a generaliza- 
tion of the response formula piI.12[) . the linear operator 
replacing Lv is no longer symmetric, in agreement with 
the breakdown of Onsager reciprocity relations. We will 
see in the next section that an inverse formulation of our 
linear response question does not suffer from the above 
complication and it also becomes remarkably easier to 
study outside the weak driving regime. 
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C. Local response to nonlocal perturbation 

In this section we discuss long range aspects in the 
inverse problem, namely, how the test potential $ that 
makes a given density p stationary, is affected by a local 
change in p. The relevance of this question and some 
other approaches are considered in Section HVl 



We start by revisiting the exactly solvable model of 
Section UlIAl The confining potential U{r) and the driv- 
ing v{r,6) — {0,f{9)/r) are always considered fixed. 
For the class of radial-angular uncorrelated densities, 
p(r, 9) = p{r) q{9), the stationary equation (jll.2[) is satis- 
fied for the test potential of the form <I>(r, 9) = ^o{r) + 
^'(0), with the radial part 



$0 = - logp - U 
and with the angular part equal to 



^ = -\0gq+ I (f-^)d^ 



(III.13) 



(III.14) 



Here 9q is an arbitrary fixed angle and the constant C is 
determined from 



C 



■fde 



(III.15) 



The latter specifies the corresponding stationary current 
field as J = (0, C/r). Away from equilibrium, i.e., 
for w — § f d9 ^ 0, the angular component "^{9) 
apparently becomes a nonlocal functional of the angular 
density q{9)^ similarly to what we have observed for the 
functional dependence p[^]. 

In the general case, we are asked to find the potential 
$ that solves the stationarity equation pi.2[) as a func- 
tion of density p for given velocity field v and confining 
potential U , say on d— dimensional space. By changing 
p — > P + Sp with J 5p{x) dx = 0, we get the linear re- 
sponse equation for G = ?7 + $ + log p, 

-V • (pVSG) = - • VSiy (III.16) 
P P 

with u = logp and J — p{v — 'VG) the stationary current. 
We recognize in the left-hand side of piI.16p the action 
of the generator 

- A/(a;)+Vlogp. V/ 

for the function / — SG. The linear operator Lp gener- 
ates a reversible diffusion in potential — logp, i.e. an equi- 
librium process. The analysis is now exactly similar as 
from (|III.12p . but with i/ replacing — V, and now restrict- 
ing to a confinement region where p{x) is approximately 
constant and maximal. The source is nonzero by the 
presence of J 7^ in the right-hand side of (|III.16p . We 



refer again to Appendix IB] for the analysis, but the con- 
clusion remains that a generic local change in density Sp 
in the confinement region requires a nonlocal adjustment 
of the potential and the corresponding response function 
derives from the Green's function for the d— dimensional 
Laplacian. 

Remark that no weak driving or small current assump- 
tion was employed in the above argument! In fact, the 
presence of the (auxiliary) reversible diffusion generated 
by Lp suggests that, even far from equilibrium, there 
are symmetries in the response functions. This is indeed 
true, see Appendix IB 31 for such a reciprocity relation. 



IV. MORE ABOUT THE INVERSE PROBLEM 

One can ask how to actually construct the test po- 
tential $ that makes a given density p stationary, i.e., 
solving pi.2p . An immediate application is found in dy- 
namical fluctuation theory; a brief review is left to Ap- 
pendix [Cl In the following we present a general proce- 
dure to solve the inverse stationary problem for a class 
of densities. Next, a variational formulation suitable for 
numerical implementation will be given. 



A. A general solution 

We are back to the set-up of Section |TT] for two- 
dimensional rotational diffusion. We restrict to those 
densities p that are everywhere bounded from zero and 
for which all (nonempty) equilevel lines, p(r, 9) — const, 
are closed curves. We also stick to trivial boundary con- 
ditions at infinity, as guaranteed by a sufficient decay of 
all the fields p, U , and v. The auxiliary velocity 

c = u-V(;7 + $ + !/), i^ = logp (IV.l) 

shares with the driving field v{r, 9) an equal vorticity, 

V X c = V X w (IV.2) 

The probability current is J = p c, and the stationarity 
condition pi.2p reads 



V • c -I- c • Vi/ = 



(IV.3) 



We only need to find the vector field c(r, 9) ; then the 
potential <I>(r, 6*) can be calculated as 



<^{x) = -v{x) - U{x) + 



(w - c) • dt (IV.4) 



7: Xq-^X 



modulo a constant, where the integral is taken along an 
arbitrary curve connecting a fixed initial point .tq with 
x. 

To determine c{r,9) solving equations (|IV.2|) ~ pV.3|) . 
we first observe that it is unique by the Helmholtz de- 
composition theorem when supplying the boundary con- 
dition that the difference c — v goes to zero at infinity. 
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Still another boundary condition has to be added in the 
case the origin is not accessible and excluded from the 
domain. 

In the following we restrict ourselves again to the two- 
dimensional plane. Equation (jIV.3|l is solved by any vec- 
tor field of the form 



B. Variational approach 

Write now the Fokkcr-Planck-Smoluchowski equa- 
tion pi.2p . considered again as the inverse stationary 
problem for the test potential $, in the form 



c(r,0)=5(Kr,0))(- 



1 dv dv 

r 89' dr 



V • ( Jo - /5V$) = 



(IV.9) 



(IV.5) 



with g an arbitrary function. The latter is fixed by con- 
dition (|Tv21), 



with Jo = Jo(p) the $— independent part of the proba- 
bility current J: 



Up) - p{v - vu) - Vp 



(IV.IO) 



V X 



(IV.6) 



that after integration over the surface enclosed by any 
equilateral curve of the density, v{r, 9) = a, and using 
Stokes' theorem yields 



Parameterizing the curve v(r, 
so that 



(IV. 7) 

a by its proper length 



d^ = 
we finally get 

.9(a) = 



rAd9) +UrJJ V rd9'dr) ' G[^i,J]^\J^i-'[J-M^i)]■[J-Mp)]dx (IV12) 



Recalling that it is an elliptic partial differential equation 
for $, its solution coincides with the minimizer of the 
quadratic functional 

J'p[i'] = ^ J pS/^ ■\/^dx + V • Jo(p)da; (IV.U) 

under the unchanged boundary conditions if present. 
This formulation is suitable for numerical computations, 
see e.g. via [35l |. 

There are other variational principles of physical im- 
portance that appear intimately related to our inverse 
stationary problem. To explain those, consider the func- 
tional 



(IV.8) 



Formulae (|IV.4|) . pV.5p . and (|IV.8|) together provide an 
explicit solution for the test potential $. 

As a check we take the example (jIII.l[) : there 

c = (0, u{r)) = V 

and the curves v = a are concentric, corresponding to the 
equipotential lines for U (r) , assuming that it is monotone 
in r. Therefore (jIV.8|) gives g{a) = —u{r)/U'{r) for a = 
-U{r) ~ log Z. 

Further, the more general example of Section UlI Al has 
c = (0, j(?')/<z(0)) with non-vanishing divergence when- 
ever q is not a constant, in contrast with the form (jIV.Sp . 
The point is that this example has a velocity field which 
is not defined at the origin, it being excluded from the 
domain. Hence, Stokes' theorem in the form (|IV.7[) can- 
not be used and a modification is needed; we omit details. 

The above solution, basically obtained by a suitable 
deformation of polar coordinates, provides a class of ex- 
amples with non- vanishing radial current. The latter is 
generally the case whenever v does not decompose into 
independent radial and angular parts; compare with the 
model of Section IIII Al 



defined for all normalized densities p and all divergence- 
less currents j, V • j = 0; see Appendix [C] for its meaning 
within the dynamical fiuctuations theory. This functional 
is manifestly positive and zero only if p and j coincide 
with the stationary density respectively the stationary 
current (for the case $ = 0). It can be used to construct 
the variational functional 



AfA = . inf g[fi,j] 

j: V-J=0 



(IV. 13) 



with minimizer equal to the stationary density. This is 
a constrained variational problem that can be solved by 
Lagrange multipliers; the solution reads 



(IV. 14) 



with $ the test potential that makes the density p sta- 
tionary, cf. (|IV.9|) . 



V- [Jo(Ai)-/iV$] =0 



(IV. 15) 



In this way, the solution to the inverse stationary prob- 
lem is an essential step in constructing the variational 
functional on densities. 

Finally, recall that equation pV.15|) is equivalent to 
the variational problem J-"^ [^'] = min with J-"^ introduced 
in (jIV.ll|l . Combining with (|IV.14ll we have the relation 



-2I[p] 



(IV.16) 
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that yields the next expression for the functional X[/i] 
(changing ^' — > 25* for convenience and integrating by 
parts): 

= sup / V* • [Jo (a*) - A*V*] dx (IV. 17) 
* J 

This unconstrained variational formula is apparently 
more useful for numerical computations than (jlV.lSP 
above. 



V. CONCLUSION 

The relation between potential and stationary density 
in mesoscopic (stochastic) systems appears to be generi- 
cally long ranged whenever there is a true nonequilibrium 
driving. That long range effect is a priori distinct from 
the long range correlations under conservative dynamics 
extensively studied before, and it occurs already for free 
particles. In models of overdamped diffusions considered 
in this paper, we have linked that long range effect to the 
slow spatial decay of the Green's function for a certain 
equilibrium diffusion process (in the first order around 
equilibrium) . 

Vice versa, a similar nonlocal change of potential is 
generically needed to create a local change in the station- 
ary density. This issue appears relevant for the inverse 
stationary problem that naturally emerges in the context 
of dynamical fluctuations and nonequilibrium variational 
principles. We have indicated how these specific issues 
become mutually related, together with comparing some 
numerically feasible schemes based on the dynamical fluc- 
tuation theory that might be of use in applications. 
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APPENDIX A: MCLENNAN THEORY OF 
STATIONARY DISTRIBUTIONS 

It is useful to see how the original McLennan reason- 
ing [24| can be used to provide some physical interpre- 
tation for the stationary density of a weakly driven dif- 
fusion. We can write formula pil.lOp in the equivalent 
form, always to linear order in e, 

r+oo 

p{x) = po{x) exp / dt {v{xt) ■ VF(a;t) - V • v{xt))l 
Jo 

(A.l) 

where (•)° denotes expectation over the equilibrium pro- 
cess (that is the process (III.1[) with v = 0) started from 



position X. (For a mathematical discussion on how to 
take the limits e ^ 0, T ^ -foo, we refer to [l^].) It 
is interesting to recognize here the linear part in the ir- 
reversible entropy flux. When the density of the test 
particles is p then the instantaneous mean work done by 
the background field v is 



V ■ [{v - W)p - Vp] dx 



as the expression between square brackets is the current 
profile at density p. Clearly, for small driving that equals 



{vVV -V ■v)pdx + 0{e^ 



In other words, to linear order in the nonequilibrium 
background, —v{x) ■ 'W{x) + V • v{x) equals the mean 
dissipated work provided the particle is at x. (Incom- 
pressibility of the background fluid can be imposed by 
letting V ■ V — 0.) That linear term is exactly what ap- 
pears in the expectation in (jA.ip . Therefore, the linear 
nonequilibrium correction to the Boltzmann distribution 
corresponds to the total dissipated work under the equi- 
librium relaxation process as started from different initial 
configurations. 



APPENDIX B: MORE TECHNICAL ASPECTS 
OF THE NONLOCALITY 



Nonlocality in (equilibrium) transient 
distributions. 



The nonlocal features as discussed in the present paper 
refer to fluctuations and responses in steady nonequilib- 
ria. Nevertheless, their origin takes us to Poisson equa- 
tions for reversible dynamics, see (jIII.12p and piI.16p . 
We therefore start here with a look at equilibrium dy- 
namics but in the transient regime. 

Take a reversible {v = 0) diffusion in a potential land- 
scape V at equilibrium, i/ =^ log/? = ~V (up to an ir- 
relevant constant). Perturb the system by changing the 
potential at time t = to V + 5V and let the system re- 
lax towards a new equilibrium. The evolved distribution 
at time t he pt = p + Spt , corresponding to the effective 
time-dependent potential v + 6i>t, 5vt = 5pt/ p. By the 
linear response theory, 

A 5vt = LvSvt + -V • ipVSV), 5vo = (B.l) 
dt p 

Using that the second term on the right-hand side equals 

1, 



P 



-V • {p\I5V) = Lv5V 



we find 



5vt^(^j e'^'^ds^LvSV 
= (e*-^^ - 1) 6V 



(B.2) 
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where the equilibrium condition has been used. As a 
result, denoting with pt{x,y) the transition kernel, 



Siytix) = ~SV{x) + / dyptix,y)SV{y) 



6vt{x) 
5V{y) 



= -^{x - y)+pt{x,y) 



(B.3) 



Remember that 6vt is only determined up to a constant, 
which explains why lim4_>oo Svt = —SV + {5V)p differs 
from the "naturally expected" value —SV; in the above 
the additive constant has been fixed by the initial condi- 
tion Si'o — 0. The nonlocal part in the linear response for 
fixed time t thus exactly equals the transition probability 
density pt{x,y). It is nonlocal in the sense that over dis- 
tances where p is approximately constant and maximal 
( or around the minimum of V ), there is only slow decay 
in \x — y\. As time grows larger, that effect typically dies 
out, restoring a strictly local response in the infinite-time 
limit. 

As an example, the standard one-dimensional 
Ornstein-Uhlenbeck process (or, oscillator process) cor- 
responding to the potential V{x) — x'^/2 has a response 
function with the large-time asymptotics 



[iyt{x) - i^t{x')] = S(x' - y) - S{x - y) 



SV{y) 



+ e-\x' - x)y p{y) -\- 0{e~'*) (B.4) 



the nonlocal component of which has a weight exponen- 
tially damped in time. 



2. Green's function in confinement region 



The response analysis in Sections IIII BfUlTC] reduces 
the problem to finding the Green's function 



Ll'^g{x,y)^-d{x-y) 



(B.5) 



with Lp — • (pV) generating diffusion in the poten- 
tial — log p. It has an explicit solution in terms of the 
transition kernel (or transition probability density)): 



G{x,y) = 



+ 00 



[Pt{x,y) - p{y)] dt 



(B.6) 



In a confinement region where p is approximately 
constant around its maximum, one expects that Q{x,y) 
is essentially determined by a free diffusion. To asses 
this conjecture, we need to understand how the inhomo- 
geneities in p outside the confinement region influence 
the Green's function inside it and, more basically, how to 
make sense to the generally ill-defined expression (jB.6p 
for the free diffusion. 



As well known [2^, for a free diffusion (related to the 
Brownian motion as \/2Bt) in dimension d, with the 
transition kernel 



p'r{^,x)^{A^t) 

7 free / 



(B.7) 



the Green's function Q^^'^'^{x) = pt{0,x) dt only exists 
in the transient case, d > 3. For d — 1,2 the divergence 
cannot be "renormalized" via a stationary density like 
in (|B.6[) as the latter does not exist. Yet, its divergent 
part is in fact a;— independent and one has in all dimen- 
sions well defined differences 



[p'r{o,x)-p'r{o,^')]dt 

{\x'\~\x\) 



2^(log|a;'|-logN) 
i7r-''/2r(d/2- l)(|.x|2-'' 



if d = 1 

if d = 2 
if d> 3 
(B.I 



Hence, the free diffusion Green's function is well defined 
up to a possibly infinite additive constant. However the 
latter becomes irrelevant due to the "dipole" character 
of the source term in piI.12p or (|III.16|) : cf. also (|B:20)) 
below. 

It remains to see in what sense the exterior of a 
confinement region enters the properties of the (true) 
Green's function. To simulate that we consider the 
(standard) diffusion in a cube [—L/2, L/2]'^ with re- 
flexive boundary conditions. The transition kernel is 
Pt{x,y) = UUi (lt{x^^y^) with 



qt{xi,yi) 



1 

T 



-y 



cos 



(Xi 



cos 



(B.9) 



For d—\ the Green's function can be obtained explicitly: 

/•OO 

G{x,y)= [pt{x,y) ~ p{y)]dt 
Jo 

{x' + y') 



10 

L 



12 



\y 



1 

2L 



(B.IO) 



Clearly, up to a correction 0{1/L) it coincides with the 
Green's function for free diffusion; moreover, the "infi- 
nite" additive constant has been regularized and fixed by 
the length of the region. 
In general one has 



Pt{0,x) 



-y 

k 



expjifc ' 



\k\h} 



(B.ll) 



with the summation over the dual lattice, 
. . . , -27r/i, 0, 2it/L, . . .. Hence, 



fe#0 



(B.12) 
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which is to be compared with the free diffusion for which, 
formaUy, 



3. "Reciprocity relations" in the inverse problem 



1 



dfc 



(B.13) 



The latter is "infrared divergent" for d < 2 and the above 
finite lattice version just provides its particular regular- 
ization (above all it provides a cut-off of the neighbor- 
hood of fc = 0). Note that an alternative (and more 
standard) way of regularizing the free Green's function 
is to add a "positive mass" : for d — 2 one obtains 



dk 



1 

2^ JO 
1 



/c|2 +£2 

uJo{u\x\) 



du 



(B.14) 



27r 



Ko{e\x\ 



with Jq and Kq the Bessel functions of the first respec- 
tively of the second kind. Its short-distance asymptotics 
is 

GeiO,x) = ~ log(e|a;|) + o(l) (B.15) 

(7 being the Euler constant). 

These calculations indicate what are the "boundary 
effects" on the Green's function in a region with approxi- 
matively constant density profile: the difference from the 
free Green's function becomes negligible on length scales 
much smaller than size of the region. Although this com- 
parison includes the removal of an infrared divergence in 
low dimensions, the latter is well "renormalizable" in the 
above sense of finite differences. Our conclusion is that 
one could for the present purposes deal just with the 
free Green's function (although as such being ill defined.) 

As a specific example we consider the Ornstein- 
Uhlenbeck process for the diffusion in a quadratic spher- 
ically symmetric potential. Its Green's function can 
be found by solving the equation (|B.5P with p{x) oc 
exp[— y(a:)] = exp(— e^|xp/2) (for simplicity we restrict 
here to a source located at the origin). Again in two 
dimension, this has a solution 



271- A/e r 



(B.16) 



up to an arbitrary additive constant. For e\x\ ^ 1, what 
we have called the confinement region, it reads 

gov{0,x)^~^log{e\x\)+o{l) (B.17) 

in agreement with either of the two above regularization 
procedures. Similarly, for dimensions d > 2 one gets 



^ou(o,.) = -^^y^^^ ^ 



,Vir) 



(B.18) 



In contrast to section ITlI B[ the analysis of section ITlI CI 
does not make any use of close-to-cquilibrium assump- 
tions while all the same, in (jIII.16[) and for the inverse 
problem, the linear response is given in terms of a re- 
versible process. A solution to (|III.16|) can be written 
as 

6G{x) =-Jdy g{x, y){y V<5z/) (y) (B.19) 



where Q{x,y) is the Green's function (jB.5|l - (|B.6p . 
Since the transition probabilities satisfy detailed bal- 
ance, the Green's function exhibits the same symmetry: 
p{x)Q{x,y) = p{y)g{y,x). 

We now propose a partial integration which is allowed 
if the current J decays sufficiently fast at infinity; using 
moreover that J is divergenceless and that 5G is physi- 
cally determined only modulo a constant, one can then 
write (|B.19I) as 



6Gix) = JdyVy 



G{x,y) 
p{y) 



[J5v){y) 



and therefore 



. „ 5G , „ 8G 



'5v{y) 



5v{x) 



(B.20) 



fB.21) 



That relation vaguely resembles an Onsager reciprocity: 
— ^^^y^ y xG can be interpreted as the extra (gradient) 
force at x needed to have a deltafunction-like response at 

y- 



APPENDIX C: DYNAMICAL FLUCTUATIONS 

Variational principles often arise in the context of fluc- 
tuation theory. For example, the entropy and related 
thermodynamic potentials play a role both as important 
equilibrium variational functionals (second law) and as 
functionals governing the equilibrium fluctuations (Ein- 
stein's fluctuation theory). Here we sketch how the func- 
tionals Q and I introduced in Section IIVBI also fit such 
a scheme; a more technical account of this problem for 
overdamped diffusions can be found in (211] . 

We consider the diffusion process (jll.ip with the test 
potential $ set to zero. For any random realization xt, 
t > 0, of this process we introduce the empirical occupa- 
tion density 



1 F 

Pt{z) = y y ^{xt - z)dt 



(C.l) 



yielding a power law decay \x\ ife|a;|<Cl 



that counts the relative time spent at each point z. Ap- 
parently, pt{z) is the random density dependent on the 
realized history. In the limit T ^ 00 it converges to the 
stationary density p, with probability one by the ergodic 
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theorem. The main result reads that for large but finite 
times T, the probability of fluctuations of pT has the 
asymptotics given by the large deviation law 

P[Pt = /i) oc e-^^I'^l (C.2) 

in which the functional of Section IIV Bl is recog- 
nized as an exponential decay rate. Then, the varia- 
tional inequality > T[p\ = is a mere consequence 
of the fluctuation law (|C.2p . As observed before, find- 
ing T[p\ amounts to solving the inverse stationary prob- 
lem (IIV.i4p - pV.i5p or, equivalently, to evaluating the 



variational expression (|IV.17p . 

Similarly, the functional (jIV.12p reveals to be the ex- 
ponential decay rate in the large deviation asymptotics 
of the joint probability law for the empirical occupation 
times and the empirical current. For proofs and for more 
details see [2l|. 

The large deviation theory for stochastic systems has 
been started and rigorously established by Donskcr and 
Varadhan, 0, H^l- In the physics literature, these meth- 
ods go back to the seminal work of Onsager and Machlup, 
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